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Abstract 


Methanol-steam reforming can be utilized as a fuel processing system for hydrogen fuel cells. A study of the reacting flow in packed-bed and 
wall-coated catalytic reactors is presented. The wall-coated reformer has a smaller power requirement for delivering fuel than the packed catalytic 
bed reformer. Also, the coated catalytic layer has a smaller thermal resistance compared to the packed catalytic bed. This yields enhanced thermal 
field management for the wall-coated reformer that is essential for reformer performance. Understanding the transport in reformers is essential for 
improving both the efficiency of the reforming process and the quality of the processed fuel. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Increased energy demand and environmental protection con- 
siderations promote the need for improved energy conversion 
systems. Among these systems, fuel cells are attractive due to 
their low pollution, high potential efficiency and ease of recharg- 
ing. Various types of fuel cells provide a variety of suitable 
solutions for different energy requirements. For the applica- 
tions of portable electrical power, especially for personal mobile 
electronics, proton exchange membrane (PEM) fuel cells are 
considered desirable owing to their relatively low operating 
temperatures, small volume and light weight [1]. In particu- 
lar, miniature PEM fuel cells based on micro scale design and 
fabrication technologies are compatible with the continuously 
growing areas of integrated circuits (IC) and microelectrocham- 
ical systems (MEMS) [2]. This yields the requirement for 
developing compatible fuel processors for converting hydrocar- 
bons to hydrogen rich gas as a fuel for miniature fuel cells. 
Methanol has been considered as a suitable hydrocarbon for 
miniature PEM fuel cell applications because of its high vol- 
umetric and gravimetric energy densities, safe handling, ease 
of storage and conversion to hydrogen at low temperatures 
(~250 °C) [2,3]. 
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The design of fuel processors or reformers can correspond to 
small volume, small weight, low cost, long lifetime, low byprod- 
uct (carbon monoxide), ease of manufacture, low flow resistance, 
low thermal resistance and short start-up time [3]. Miniature or 
micro methanol reformers provide the potential to achieve these 
requirements and have attracted great interest [4]. Several micro- 
reformers based on MEMS technologies have been fabricated 
and tested [5—14]. Micro-reformers made by packing catalyst 
particles into micro channels have the advantage of easy cata- 
lyst loading. However, high flow resistance resulting from the 
packing catalyst particles may limit the design and performance 
of micro-reformers [15-17]. 

Recently, micro-reformers with wall-coated or suspended 
catalytic layer configurations have been studied due to their 
lower transport resistances compared to packed-bed reformers 
[5-11]. In the present work, studies of the reacting flow in a 
packed-bed and a wall-coated reformers are carried out. 


2. Experiment 
2.1. Packed-bed reformer 

In previous studies [17,18], experiments were carried out in a 
packed-bed reformer with the BASF K3-110 CuO/ZnO/A1203 


based catalyst. Catalyst pellets were ground into powder and 
sieved to obtain an average diameter of 25 microns. Borosilicate 
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Nomenclature 

A lst pre-exponential coefficient [m3 (s kg)~!] 

B 2nd pre-exponential coefficient [m3 (s kg)~!] 

Cb molar concentration of the stream 
[kmol m~3] = p/(RT) for ideal gas 

E activation energy [kJ kmol~!] 

Fm methanol inlet flow rate [mol s7 !] 

kp forward reaction rate constant of methanol 
decomposition reaction [kmol (m? s)~!] 

kR forward reaction rate constant of steam-methanol 
reforming reaction [ls~ 1] 

Meat mass of catalyst [kg] 

M molecular weight [kg kmol~!] 

p pressure 

i" rate of production [kmol (m? s)~"] 

R universal gas constant 

SMR steam to methanol molar ratio 

T temperature [K] 

u axial direction gas mixture velocity [m s7 !] 

Vin volumetric flow rate of methanol and water mix- 


ture at the inlet of the reformer 
density of catalyst [kg m7°] 
Xi mole fraction of species i 


Greek letter 


p density of gas mixture [kg m~7] 

E reaction progress parameter 

Subscripts 

1 methanol 

2 water 

3 hydrogen 

4 carbon dioxide 

5 carbon monoxide 

b bulk value 

D value for methanol decomposition reaction 
R value for steam-methanol reforming reaction 


glass tubes with 1.5 mm inner and 1.8 mm outer diameters were 
filled with the catalyst powder and the two ends of the pack- 
ing bed were held by inserted deactivated silica wools. Fig. 1 
shows a schematic of the experimental apparatus. Mixtures of 
distilled water and methanol with molar ratio 1.1 flowed into the 
evaporator at specified flow rates using a syringe pump (Genie 
programmable syringe pump, Kent Scientific Corp.). A k-type 
thermocouple was inserted in the middle of the catalyst packing 
bed to obtain a characteristic reformer temperature. The reac- 
tants were heated and the gas then flowed into the reformer. The 
reformer was heated to supply the energy for the endothermic 
reforming reaction. The temperature of the outer surface of the 
reformer was controlled by utilizing a temperature controller 
(Minco CT-16A). The gas leaving the reformer passed through 
a cold trap to separate the unreacted reactants, i.e. methanol 
and water. The total flow rate of the product gas mixture (H2, 
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Fig. 1. Apparatus of experimental setup. 
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CO and CO) at the outlet of the cold trap was measured with 
a flowmeter (ADM 2000, Agilent). Tests were carried out first 
for a specific inlet flow rate of methanol and water mixture and 
the temperature was then varied for different runs [19]. The 
sequence was then repeated for a new inlet flow rate. 

The concentration of H2 of the product gas mixture was 
measured with a H? sensor (Robust Hydrogen Sensor, DCH 
Technology Inc.) and the flow rate of H2 was then obtained by 
multiplying the total product gas mixture flow rate and the H2 
concentration. Assuming a small CO concentration, the flow rate 
of CO7 is then obtained by subtracting the H2 flow rate from the 
total flow rate. 

Amphlett et al. [20] carried out experiments in a packed-bed 
reformer with methanol and water mixtures with the BASF K3- 
110 catalyst. They obtained the constants in the Arrhenius law 
expressions which provide the semi-empirical relations for the 
methanol reforming and decomposition reaction rate constants, 
kp and kp, respectively. Utilizing the reaction rate constants of 
Amphlett et al. [20] in conjunction with a “reaction progress 
variable” (cf. Appendix A), Park et al. [17,18] obtained theo- 
retical results for the methanol conversion, the product mixture 
flow rate and the molar concentrations of CH30H, H20, Ho, 
COz and CO in a packed-bed steam-methanol reformer at atmo- 
spheric pressure. Comparisons between predicted and measured 
Hə and CO3 flow rates were made and the results were in good 
agreement. 

In the present work, another CuO/ZnO/A1203 based catalyst 
(BASF F3-01) was used in both a packed-bed and a wall-coated 
reformer. Tests of the catalyst were first carried out in a packed- 
bed reformer. Specifically, flow rates and concentrations of H2, 
CO, and CO were measured at various reacting temperatures 
and flow rates of the methanol and water mixture. Particle size 
is an important factor for the catalyst because of intraparticle 
diffusion [28]. Therefore, two different catalyst particle sizes, 
75 and 150 microns, were tested. 15 to 16mg of the BASF 
F3-01 catalyst particles were loaded into borosilicate glass 
tubing having an inner/outer diameter of 1.5 mm/1.8 mm, 
respectively. The loading length of the catalyst was about 
1.1 cm of the 75 micron catalyst particles and about 0.9 cm for 
the 150 micron particles. The two ends of the catalyst packing 
bed were sealed with silica wool. A temperature controller was 
used to maintain the outer surface temperature of the tubing at 
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Table 1 

Reaction progress variable (1 =CH30OH, 2 = H20, 3 =H2, 4=COz, b= bulk) 

é Xib X2b X3b X4b Up My 

0 1 SMR 0 0 ha 18SMR + 32 
SMR + 1 SMR + 1 SMR + 1 

(2) 1-é SMR —& 3E E (SMR + 1 + 2Ẹ)up,in 18SMR + 32 
SMR + 1 +24 SMR + 1 +24 SMR + 1 + 2& SMR + 1 +2 SMR + 1 SMR + 1 + 2& 

i 0 SMR- 1 3 1 (SMR + 3)ub,in 18SMR + 32 

SMR + 3 SMR + 3 SMR + 3 SMR + 1 SMR + 3 


the middle of the catalyst packing bed. Heating upstream of the 
reformer provided a gas phase flow of methanol/water mixtures 
at the inlet to the reformer. The experimental procedure was the 
same as described previously. The detailed operating conditions 
of these tests are listed in Table 2. 

Deactivation of the CuO/ZnO/Al,03 based catalyst in 
packed-bed reformers, especially the degradation caused by 
sintering, is expected even in the recommended working tem- 
perature range [29,30]. Therefore, it is essential to test the 
repeatability of the data before carrying out the kinetic tests. The 
repeatability tests for BASF F3-01 were carried out at 230°C 
and a 20 pL min! inlet flow rate of a water and methanol mix- 
ture with the steam to methanol molar ratio SMR=1.1. Fig. 2 
shows the repeatability test results for the conversion which is 
defined as 


; X1b,in — *1b,out 
conversion = ———————— 


x 100 (1) 
X1b,in 

The subscript 1 refers to methanol and the subscript b refers to 

the bulk value. For the determination of the mole fractions x1p in 

and X1p,out refer to Appendices A and B and Table 1. 

The amplitude of the reduction rate of the conversion, which 
corresponds to the degradation rate of the catalyst activity, 
decreases gradually with respect to testing time for both cat- 
alyst particle sizes. This result is in agreement with the trend 
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Fig. 2. Results of repeatability test for packed-bed reformers, testing at 230°C, 
inlet flow rate of 20 uL min—!, SMR = 1.1, catalyst loading = 16 mg. 


of deactivation of the CuO/ZnO/Al203 based catalyst reported 
previously [29,30]. After 4h of reforming, the conversion was 
reduced from 43 to 38%. All other tests in this work were 
done within 4h for each reformer newly loaded with catalyst. 
It is noted that Thurgood et al. [30] reported that the rate of 
decline of the catalyst activity increases with increasing reac- 
tor temperature. Results shown in this work were obtained 
from the experiments performed at temperatures not higher 
than 230°C, 

Tests were also carried out to obtain the CO concentrations 
of the product gas mixtures. The experimental conditions and 
procedures were the same as described in the previous section. 
With a specified inlet flow rate of the methanol-water mixture 
and reformer temperature, the concentration of H2 was mea- 
sured with the DCH H3 sensor and the product gas mixture was 
collected with a gas sampling bag (SKC Inc.) after the prod- 
uct gas mixture flow rate was stabilized. The collection of the 
product gas mixture for each test was completed within 3 h. The 
product gas mixture flow rate was checked every hour during the 
test to check the stability of the system. The CO concentrations 
in the collected gas sample bags were measured by utilizing 
a Fourier Transform Infrared Spectrometer (Magna 760 FTIR, 
Nicolet Instrument Corp.). 


2.2. Wall-coated reformer 


Coating of the catalytic layer was carried out following the 
wash-coated procedure developed by Bravo et al. [31]. Borosili- 
cate tubes with the same size as used for the packed-bed reformer 
were wash-coated with the slurries made by mixing and dispers- 
ing the BASF F3-01 catalyst and distilled water. The length of the 
catalytic coating layer varies from 2 to 6cm. The average coat- 
ing thickness was about 100 wm. This value was obtained from 
the scanning electronic microscopic (SEM) images of the wash- 
coated reformers as shown in Fig. 3. The average loading of the 
catalyst was about 8 mg. Experiments and measurements were 
carried out with the same testing sequence as for the packed-bed 
reformer tests. Outer wall surface temperatures at the locations 
of 1/4, 1/2 and 3/4 of the coating length were monitored and 
maintained at the desired temperature. The operating conditions 
were listed in Table 2. Fig. 4 shows the result of the repeatabil- 
ity test of the wall-coated reformer. After 4h of reforming, the 
conversion was reduced from 28% to 20%. All other tests in this 
study were done within 4h for each reformer newly loaded with 
catalyst. 
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Fig. 3. Scanning electronic microscope images of wall-coated reformer. Tube i.d. = 1.5 mm. Catalyst coating thickness is about 100 wm. 


Table 2 

Reformer test conditions 

Water to methanol molar ratio (SMR) 1.1 

Feeding rate 2-30 pL min”! 
Reacting temperature 200-250 °C 


Environment pressure l atm 
Environment temperature Room temperature 


3. Results and discussion 
3.1. Methanol conversion 


Comparisons of the conversion of packed-bed and wall- 
coated reformers at three different reforming temperatures are 
shown in Fig. 5(a)—(c). The low conversion rates are mainly due 
to the low reactor temperatures and the small amount of cata- 
lyst loading; the methanol conversion increases with increasing 
reactor temperature or increasing the amount of catalyst loading 
(Meat). It is shown that the conversion of the packed-bed reform- 
ers differs only slightly for the two different particle sizes (75 
and 150 um). Purnama et al. [28] carried out experiments on 
two packed-bed reformers with a CuO/ZnO/A1203 Siid-Chemie 
catalyst for particle sizes from 0.71 to 1.0mm and from 0.45 to 
0.5 mm, respectively, and showed that the smaller catalyst par- 
ticles produced higher conversion. Jiang et al. [23] carried out 
methanol reforming experiments on CuO/ZnO/Al,03 BASF S3- 
85 catalyst particles for sizes from 150 to 590 microns and noted 
no effect of particle size on methanol conversion. The catalyst 
particle sizes utilized in the present study were much smaller 
than were used by Purnama et al. [28] and were comparable to 
those used by Jiang et al. [23]. Consistent with the results of 
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Fig. 4. Results of repeatability test for wall-coated reformer, testing at 230°C, 
inlet flow rate of 20 pL min7!, SMR=1.1. 


Jiang et al. [23], there was found to be little effect of the catalyst 
particle size on the conversion. 

Bravo et al. [31] carried out methanol reforming experiments 
on (a) packed-bed and (b) wall-coated reformers in 4.1mm 
diameter tubes that were (a) packed with 100 to 250 micron cata- 
lyst particles or (b) coated with catalyst on the inner wall surface 
of the tubes. They used CuO/ZnO/A1203 BASF F13456 catalyst. 
At the reforming temperature of 230 °C, their results showed that 
the wall-coated reformer produced much greater conversion than 
that of the packed-bed reformer. In the present study (with the 
catalyst BASF F3-01), the conversions produced by wall-coated 
reformers are very close to those of the packed-bed reformers. 


conversion (%) 


(b) 


conversion (%) 
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conversion (%) 


Fig. 5. Comparison of conversions of packed-bed and wall-coated reformers. 
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Karim et al. [16,32] carried out experiments and transport 
analyses on packed-bed and wall-coated reformers. Their exper- 
imental results on packed-bed reformers with tube diameters 
1, 1.75 and 4.1 mm show that the conversion increased with 
decreasing diameter [32]. In their analysis, they noted a signifi- 
cant temperature difference in the packed-bed reformer with the 
temperature near the core being much lower than the nominal 
reforming temperature (which is measured at the wall surface); 
the reactions are endothermic and the heat was supplied to the 
reformer from the outside surface of the wall [16]. This results in 
the average temperature of the catalyst bed being significantly 
lower than the temperature near the wall and the temperature 
difference between the surface and the core of the packed-bed 
reformer increased with increasing tube diameter. The exper- 
imental results of wall-coated reformers with tube diameters 
0.53, 1.75 and 4.1 mm obtained by Karim et al. [32] showed 
that the conversions are almost the same. Their analysis of the 
wall-coated reformer with a tube diameter of 0.53 mm showed 
that the temperature difference in the thin catalyst coating layer 
is very small so that the average temperature of the wall-coated 
catalyst layer is very close to the wall temperature. Additional 
differences between the packed-bed and wall-coated reformers 
should also be present corresponding to the different mass trans- 
port and surface reaction characteristics. Experimental results of 
Karim et al. [32] showed that at the same wall temperature the 
conversion produced from the packed-bed reformer with tube 
diameter 1 mm is greater than from the wall-coated reformers 
with tube diameters 0.53, 1.75 and 4.1 mm. They postulated this 
to result from the catalyst structure and properties that changed 
during the coating process. 

It would appear that if the conversion is dominated by the 
temperature effect (distribution) then for the same wall tem- 
perature the higher average catalyst layer temperature of the 
wall-coated reformer should yield greater conversion than that 
from the packed-bed reformer. This is apparently the condition 
for Bravo et al. [31] but not for the present experiments which 
show only little difference in conversion between the packed- 
bed and wall-coated reformers at the same wall temperature. In 
the present study, a 1.5 mm diameter tube was utilized for both 
packed-bed and wall-coated reformers and it is expected that for 
the packed-bed reformer that the temperature difference in this 
1.5 mm tube is less than that for the larger 4.1 mm diameter tube 
that was used by Bravo et al. [31]. Thus, in the present experi- 
ments the average temperatures of the packed-bed reformer and 
wall-coated reformer should be closer to each other than for 
the larger diameter tube of Bravo et al. [31] and the effect of 
the temperature differences, although present, should be dimin- 
ished. It is noted that the experimental results obtained by Karim 
et al. [32] showed the conversions of packed-bed and wall-coated 
reformers with a tube diameter of 1.75 mm to be almost the same 
at the same wall temperature. 


3.2. Carbon monoxide 
Fig. 6 shows the results for the carbon monoxide con- 


centration produced in the packed-bed reformers filled with 
150 micron catalyst particles. The concentration of CO increases 
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Fig. 6. Comparison of carbon monoxide concentrations produced from packed- 
bed reformer at various reforming temperatures. 


with respect to the value of meat/Vin which is proportional to 
the residence time or contact time that the reactants and prod- 
ucts are in the reformer. At the same value of Mmeat/ Vin, the CO 
concentration is greater at higher temperature. Park et al. [33] 
measured the CO concentration produced in a micro packed- 
bed reformer with a specified catalyst loading weight and inlet 
methanol-water mixture flow rate for different reforming tem- 
peratures. They also showed an increasing concentration of CO 
with increased temperature. Agrell et al. [21] reported a decreas- 
ing concentration of CO with increasing steam to methanol ratio 
(SMR). The benefit of utilizing a greater value of SMR to mini- 
mize CO generation is offset by the lower hydrogen generation 
rate that results because methanol contributes more hydrogen 
per mole than water. Accordingly, a compromise value of SMR 
(1.1) was utilized in the present system. 

Amphlett et al. [20] utilized the methanol decomposition 
reaction as the model to predict the concentration of CO and 
obtained a result for the CO concentration that is linearly 
dependent on the residence time. Recent studies [21,27,34] 
have suggested that in steam-methanol reforming with the 
CuO/ZnO/A1203 based catalyst, CO is produced as a secondary 
product from the reverse water-gas-shift reaction (Eq. (2)) 


CO2 +H» — CO + HO (2) 


rather than from methanol decomposition. 

Based on their experimental results, Purnama et al. [28] fur- 
ther suggested that the controlling kinetics of CO formation is 
dependent on the magnitude of the residence time and reform- 
ing temperature. The experimental results of CO concentration 
obtained in the present study are consistent with the trend pre- 
sented in the study by Purnama et al. [28]. 

Fig. 7 shows a comparison of the CO concentration produced 
from two packed-bed (filled with catalyst particles of diameters 
150 microns and 75 microns) and wall-coated reformers at the 
reforming temperature of 230°C. It is shown that the CO con- 
centrations produced from all the tests are similar to each other 
at the same value of MMeat/Vin. 
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Fig. 7. Comparison of carbon monoxide concentrations produced from packed- 
bed and wall-coated reformers. 


Bravo et al. [31] carried out similar tests in a wall-coated 
reformer with a 25 micron thick catalyst layer wash-coated on 
the inner wall of a 4.1 mm inner diameter glass tube. They also 
carried out tests in a packed-bed reformer with a 4.1 mm i.d. 
glass tube filled with 100-250 micron catalyst particles. The CO 
concentration was higher for the wall-coated reformer than for 
the packed-bed reformer at the same value of the catalyst weight 
to methanol inlet flow rate ratio (Meat/Fm). It is noted that Bravo 
et al. [31] presented their CO results with respect to the conver- 
sion. Using this result, and their conversion results with respect 
tO Meat/ Fm, we obtained the result noted above for CO concen- 
trations with respect to mcat/Fim. It is also pointed out that Bravo 
et al. [31] showed that the CO concentration is less for the wall- 
coated reformer than for the packed-bed reformer at the same 
value of the conversion. These differences of CO concentra- 
tions between the packed-bed and wall-coated reformers were 
attributed to result mainly from the temperature variation across 
the reformer as described in the previous section. 

In the present study, however, the differences of the CO 
concentrations for the two reformers are not as large as were 
observed by Bravo et al. [31]. We have pointed out that the tube 
being used as the reformer in the present study is smaller than 
that used by Bravo et al. [31]. The temperature difference across 
the catalyst bed in the packed-bed reformer should be smaller in 
the smaller tube. Thus, the average temperatures of the packed- 
bed and wall-coated reformers are closer to one another and the 
effect of the temperature differences is diminished in the present 
study. 


4. Conclusions 


Packed-bed and wall-coated methanol-steam reformers were 
studied in the present work. The results show that the methanol 
conversion was reduced by 5-8% after 4 h of continuous reform- 
ing and the amplitude of the reduction rate decreased gradually 
with testing time. Under the same reforming condition, the pack- 
bed reformer filled with 150 micron or with 75 micron catalyst 
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particles produced the same methanol conversion and CO con- 
centration. At the same reformer wall temperature, the methanol 
conversion and CO concentration produced from the packed- 
bed reformer and from the wall-coated reformer are the same 
for the same value of the catalyst mass to inlet fuel flow rate 
ratio (Meat/ Vin). This result indicates that at the same wall tem- 
perature that the average temperature of the catalyst bed in the 
packed-bed reformer is close to the average temperature of the 
catalyst layer in the wall-coated reformer. 
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Appendix A. Reaction progress variable (Park et al. 
[17,18]) 


The kinetics of steam-methanol reforming reaction has been 
studied intensively; a number of kinetic models have been pro- 
posed [20-27], and the model of Amphlett et al. [20] is utilized. 
In Amphlett’s analysis, the steam-methanol reforming reaction 
can be modeled with two major reactions: 


Reforming reaction CH30H + H20<=>3H} +CO: (Al) 
R’ 


Decomposition reaction CH30H- or, + CO (A2) 


The water-gas shift reaction was assumed in Amphlett’s anal- 
ysis to have negligible impact. Reaction rates of the reforming 
and decomposition reactions can be expressed as follows: 


ip = keep (A3) 
m ER 
kg = we [Ar + Br In(SMR)] exp Pag (A4) 
iD = kp (A5) 
Ep 
EE EEA |- z (A6) 


The reaction rate of the reforming reaction is found to be 
much greater than the decomposition rate due to the selectivity of 
the catalyst and is thus the main source of hydrogen generation. 
To obtain the amount of carbon monoxide generation, Amphlett 
used the decomposition reaction and is utilized in this study. 

If initially there is one mole methanol and SMR moles of 
water vapor, the molar fractions of methanol and water vapor 
can be expressed as x1p = 1/(1 + SMR) and x2 = SMR/(1 + SMR), 


respectively. There is no hydrogen or carbon dioxide at the 
entrance to the reformer. 

As reaction occurs on the catalyst as reactants flow toward the 
downstream of the reformer, hydrogen is generated as the main 
product and carbon dioxide is produced as the main byprod- 
uct. Presuming €, moles of methanol having been reacted after 
some distance z from the reformer entry, referring to the steam- 
methanol reforming reaction (A1), £, moles of water vapor have 
also been reacted and 3€, moles of hydrogen as well as £, moles 
of carbon dioxide have been generated. The molar ratio of each 
of the four major species can then be obtained as a function of £, 
and is listed in Table 1. The value of £, increases as the reforming 
reaction continues. With methanol being completely reformed, 
there will be 3 moles of H? and 1 mole of CO% as products; 
SMR-1 moles of water and no methanol are left as reactants, as 
listed in Table 1. 

With knowing the molar fractions of the four major species, 
the molecular weight, Mp, of the reacting gases flow as function 
of £, can therefore be obtained. Also, by assuming ideal gas, 
isothermal and constant pressure, the expressions for average 
velocity, up can be derived, as summarized in Table 1. 

The bulk density, pp, can also be obtained from the ideal gas 
law: 


b 
Pb = Mp (=) = Mbc (A7) 


Appendix B. Procedure for evaluating mole fractions 


Refering to Appendix A and Table 1 the mole fraction x1» is 
given by 


1=£ l= 
Xib = 


—SME+142 314-22 


(Bl) 


where SMR = 1.1. To obtain xjp, out, the value of Eout is obtained 
from Table 1 for up = Up out: (B2)up, out = ae LY The mea- 
sured values of Up in and Up out then yield the value of Eout Which 
then yields x1b,out from Eq. (B1). 


Note that (B3) 


Xibin = 57 


1 
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